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Abstract, interacting Dark Energy models have been introduced as a possible alternative to the 
standard ACDM concordance cosmological scenario in order to ease the fine-tuning problems of the 
cosmological constant. However, the interaction of the Dark Energy field with other massive parti- 
cles in the universe induces also an effective modification of structure formation processes, leading 
to a different dynamical behavior of density perturbations with respect to the standard scenario. In 
particular, high-resolution N-body simulations have recently shown that also the structural properties 
of highly nonlinear objects, as e.g. their average concentration at a given mass, could be signifi- 
cantly modified in the presence of an interaction between Dark Energy and Dark Matter. While a 
constant interaction strength leads to less concentrated density profiles, a steep growth in time of the 
coupling function has been shown to determine a large increase of halo concentrations over a wide 
range of masses, including the typical halos hosting luminous spiral galaxies. This determines a sub- 
stantial worsening of the "cusp-core" tension arising in the standard ACDM model and provides a 
direct way to constrain the form of the Dark Energy interaction. In the present paper we make use 
of the outcomes of some high-resolution N-body simulations of a specific class of interacting Dark 
Energy models to compare the predicted rotation curves of luminous spiral galaxies forming in these 
cosmologies against real observational data. Our results show how some specific interacting Dark 
Energy scenarios featuring a steep growth in time of the coupling function - which are virtually in- 
distinguishable from ACDM in the background - cannot fit the observed rotation curves of luminous 
spiral galaxies and can therefore be ruled out only on the basis of dynamical properties of small-scale 
structures. Our study is a pilot investigation of the effects of a Dark Energy interaction at small scales, 
and demonstrates how the dynamical properties of visible galaxies can in some cases provide direct 
constraints on the nature of Dark Energy. 
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1 Introduction 

The presently accepted standard cosmological model based on a cosmological constant A and on 
Cold Dark Matter (CDM) particles as the two dominant components of the Universe has proved to 
be extremely successful in describing the overall evolution and the structural properties of the Cos- 
mos on large scales. Ever improving datasets have allowed in the last decade to tightly constrain the 
properties of Dark Energy (DE) at low redshifts and to show that the observed accelerated expansion 
of the Universe [1-3] is driven by some field with an equation of state extremely close to wde = — 1 
at the present epoch [see e.g. 4-7]. Furthermore, the success of the CDM paradigm in explaining the 
formation and the statistical properties of the large-scale structures observed in the Universe has been 
strongly supported by the parallel improvement of observational surveys and numerical simulations. 
Nevertheless, the extreme fine-tuning of the cosmological constant value required to fit the data on 
the background evolution of the Universe poses a serious naturalness problem to the model [8] and 
represents the main motivation for the investigation of alternative scenarios. Additionally, a signifi- 
cant number of astrophysical observations at small scales seem to show discrepancies with respect to 
the predictions of the ACDM model. These range from the low number of detected luminous satel- 
lites in galactic-size CDM halos (the "satellite problem" [9, 10], but see also e.g. [11] for a possible 
explanation), to the observed low baryon fraction in galaxy clusters [12, 13], the large peculiar ve- 
locities detected in the bulk motion of galaxies [14] or in systems of colliding galaxy clusters as the 
"Bullet Cluster" [15], to the shallow observed density profiles of CDM halos (the so-called "cusp- 
core" problem) of dwarf galaxies [16, 17], spiral galaxies [18, 19], and galaxy clusters [20, 21]. In 
the present paper, we will consider the latter problem in relation to the observed dynamical properties 
of spiral galaxies and discuss the impact that some specific alternative cosmological scenarios might 
have on this issue. In fact, based on the shortcomings of the ACDM cosmology concerning both its 
fundamental nature and its observational properties at small scales, a significant number of alterna- 
tive models have been proposed in recent years, ranging from DE scenarios based on the dynamical 
evolution of a classical scalar field as for the case of quintessence [22-24] or k-essence [25], to spe- 
cific modifications of General Relativity at cosmological scales [26], to Warm Dark Matter (WDM) 
scenarios [27, 28]. A particularly interesting class of alternative models is given by interacting DE 
scenarios [29, 30], where a direct exchange of energy-momentum between the DE field and massive 
particles takes place. Such models might provide early scaling solutions for the DE density thereby 
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alleviating the fine-tuning issues of the cosmological constant. Additionally, they predict the exis- 
tence of a fifth-force between massive particles coupled to the DE field, which determines significant 
and potentially observable effects on the growth of cosmic structures [see e.g. 31-33] and is there- 
fore likely to have a direct impact on the small scale failures of the ACDM scenario. However, the 
phenomenology of coupled DE models (cDE) has been shown to be extremely diverse depending 
on the nature of the interaction [34] and while providing a fully viable alternative to ACDM for a 
wide range of realizations, it might in some specific cases also significantly worsen the problems of 
the ACDM scenario on small scales while being still in full agreement with background observables. 
This provides a direct way to strongly constrain the parameter space of cDE models and represents 
the main focus of the present work. More specifically, we aim at quantifying the impact that a steeply 
growing coupling between DE and CDM particles has on the predicted rotation curves of luminous 
spiral galaxies based on an NEW density profile, due to the change induced by the interaction on 
the normalization of the Concentration-Mass relation of their host DM halos. The latter effect has 
been investigated through numerical N-body simulations of structure formation both for the case of 
a constant [32] and a time-dependent [34] interaction strength, and we seek here a direct comparison 
of these results with observational data from galaxy rotation curves. This will allow us to rule out 
some specific realizations of the cDE scenario based on rotation curves data only. 

The paper is organized as follows. In Section 2 we summarize the main features of cDE models 
and we highlight the specific scenarios considered in the present work. In Section 3 we discuss how 
we model the mass distribution within spiral galaxies and how the velocity profiles are derived based 
on the structural properties of the host CDM halo. In Section 4 we provide constraints on the cDE 
models under investigation from our sample of galaxy rotation curves, while in Section 5 we discuss 
our findings and draw our conclusions. 

2 Coupled Dark Energy models and Simulations 

Coupled DE models are based on the dynamical evolution of a classical scalar field whose energy 
density = (jP' /2 + V{<j)) evolves in time according to the dynamic equation: 



that includes on the right-hand side a source term for the exchange of energy-momentum with a 
different cosmic fluid of density pc- In Eq. 2.1, an overdot indicates a derivative with respect to the 
cosmic time t, H = a/a is the Hubble function, Mp\ = 1/ VSirG is the reduced Planck Mass, V{(p) 
is the scalar field self-interaction potential, and the coupling function Pcifp) defines the strength of the 
DE interaction. For the purpose of the present paper, the fluid coupled with the DE scalar field will 
be always assumed to be represented by CDM particles, although cDE models with interactions to 
massive neutrinos have also been proposed as a possible solution to the cosmic coincidence problem 
[see e.g. 35, 36]. In order to preserve General Covariance, Eq. 2.1 requires the presence of an 
analogous interaction term also in the CDM continuity equation, which reads: 



while the remaining continuity equations for the baryonic and relativistic components of the universe 
are unaffected as the coupling involves only DE and CDM: 




(2.1) 




(2.2) 



Pb + 3Hpb = , 

Pr + 4:Hpr = . 



(2.3) 
(2.4) 
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As a consequence of the interaction, and assuming the conservation of the CDM particle number, 
Eq. 2.2 implies a time evolution of the CDM particle mass due to the dynamic nature of the DE scalar 
field, according to the equation: 



me(a) = mc(ao)e^v273//3c(<^)#/Mp, ^ ^2.5) 

where qq is the cosmic scale factor at the present time. The cosmological dynamics of the system 
allows for scaling solutions during matter domination that feature an Early Dark Energy component, 
thereby alleviating the fine-tuning problems of the cosmological constant. Such scaling solutions have 
been identified analytically for the case of constant coupling functions /3c((A) = Pc = const., and 
numerically for the case of variable couplings Pc{4') considered in the present work , and are known 
as the "(p-Matter Dominated Epoch" [or i;^MDE, see 30] and the "GTOv/mg-cl)-Matter Dominated 
Epoch" [or G-(/)MDE, see 34], respectively. However, if the coupling exhibits a very steep growth in 
time (as for the case of the specific models discussed in this work) the early DE scaling is absent and 
the model is therefore indistinguishable from ACDM at the background level [see again 34]. 

At the linear perturbations level, the growth of density perturbations is modified by the inter- 
action between DE and CDM through the appearance of two new terms that directly depend on the 
DE-CDM coupling in the perturbed linear- equations. The dynamic equation for linear CDM den- 
sity perturbations in the Newtonian gauge and in Fourier space is in fact modified from its standard 
form and reads [see e.g. 30-32, 37, for a full derivation of the linear perturbation equations in cDE 
models] : 
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Pbh + PcSc (1 + ^/3^^(0) 



(2.6) 



where the overdensity (5h ^ of baryons and DM is defined as ^ = ^Pb,c/Pb,c ■ The additional contri- 
bution appearing in the first term on the right hand side of Eq. (2.6) is an extra friction associated with 
momentum conservation in cDE models [see e.g. 38, for a discussion on the effects of the friction 
term] while the additional term in the second squared bracket includes the effect of the fifth-force 
mediated by the DE scalar- field for CDM density perturbations. Qualitatively, while the fifth-force 
term is always attractive and accelerates the growth of density perturbations for any value of the cou- 
pling and for any dynamic evolution of the DE scalar field 0, the friction term enhances the growth 
only when /3c(^ > 0, while it slows down the evolution of perturbations when pcfp < 0. Although in 
the present work we will be interested in specific cDE models for which the former condition always 
holds, a change in sing of the friction term can have particularly interesting effects on structure for- 
mation processes, especially at high redshifts [39, 40]. Furthermore, even for the case of cDE models 
with always positive values of I3c{(t>)4>, the friction term has been shown to determine a suppression 
of the growth of density perturbations once these become nonlinear", due to the presence of non-radial 
velocities [32, 38]. 

In the present work, we will focus on cDE models with an exponential self-interaction potential 
and with variable couplings in the form of a power of the cosmic scale factor a: 

/3e(<A(a)) = Poa^' (2.7) 

that have been proposed and studied by [34]. In particular, we will consider two of the models dis- 
cussed in [34], featuring a rapidly growing coupling function with (/3o, = (0.5, 2) and (0.75, 3), 
named "EXP010a2" and "EXP015a3", respectively. Such models have been shown to be fully consis- 
tent with present cosmological data on the background expansion history and cluster number counts. 
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Figure 1. The ratio of the halo concentration in cDE models with respect to the standard ACDM scenario as a 
function of the halo virial mass M200, as obtained from the N-body simulations of [34]. The curves represent 
the Mean (solid) and the Median (dashed) concentration within each of the four logarithmic mass bins in which 
the sample has been subdivided. 

while determining a significant increase of the small-scale matter power only at very low redshifts. 
In the present paper, we are seeking for observational constraints on this specific class of models 
based only on the observation of rotation curves of luminous spiral galaxies, without invoking any 
additional observational data. 

We will base our analysis on the numerical results of [34], that made use of the specific mod- 
ification of [32] of the widely used parallel N-body code GADGET [41] to run high-resolution hy- 
drodynamical N-body simulations of a series of cDE models with variable coupling, including the 
EXP010a2 and EXP015a3 models under investigation here. Such modified numerical code includes 
all the specific effects related to the DE-CDM interaction that have been briefly discussed above, and 
allows to follow the evolution of a cosmological volume of the universe in the context of different 
models with a common normalization of the basic cosmological parameters, and to compare the in- 
dividual and statistical properties of collapsed objects fomiing in each scenario. In particular-, we will 
consider here the results of hydrodynamical N-body simulations within a periodic cosmological box 
of 80 comoving Mpc/h aside filled with 512^ CDM and 512^ gas particles, for a mass resolution of 
2.41 X 10*^ M0//1 and 4.82 x 10'' Mq/H, respectively, and a spatial resolution of 3.5 kpc/h. By com- 
paring the outcomes of such simulations, [34] found that - differently from standard cDE scenarios 
with constant coupling - variable coupling models might determine in some specific realizations a 
significant increase of the normalization of the "concentration vs. Mass" (c-M) relation for collapsed 
virialized objects (see Fig. 1). In particular, this happens for models where the coupling steeply grows 
at low redshifts as it is the case for couplings scaling as a positive power of the scale factor, due to 
the absence in such models of a G-(/)-MDE phase [see 34, for a detailed discussion of this effect]. 
An overall increase of the normalization in the c-M relation might have a significant impact on the 
predicted rotation curves for luminous spiral galaxies and provide a powerful way to strongly con- 
strain the parameter space for this class of cDE models. More specifically, an increase by a factor 
of ~ 1.85 [as predicted by 34, for the model EXP015a3] of the expected average concentration of 
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CDM halos in the mass range of a few x 10^^ MQ/h will clearly exacerbate the marginal tension of 
observed rotation curves with the predicted radial mass distribution of a Navarro-Frenk- White [NFW 
9] density profile in ACDM, known as the "cusp-core" problem, thereby providing a direct way to 
constrain the model using only galaxy rotation curves data and without resorting on any additional 
observational dataset. 

In this paper we will adopt the outcomes of [34], where simulations were normalized in order 
to have the same amplitude of linear density perturbations at z = 0, i.e. all cosmological models were 
normalized to the same value of a^{z = 0). Recently, a new set of very large N-body simulations 
for cDE models [the CoDECS project, see 42] has been released, featuring a different choice for the 
power spectrum normalization with a common amplitude of density perturbations at the redshift of 
CMB, zcMB ~ 1100. The choice to consider the numerical results of [34] instead of [42] for our 
analysis should be considered "conservative" as a high-z normahzation of density perturbations is 
expected to result in an even stronger increase of halo concentration at a fixed mass for the same 
cosmological model. 



3 Mass modelling the Spiral's Rotation Curves 

In spiral galaxies the presence of a large amount of unseen matter, distributed very differently from 
the stellar- and HI gas disks, is well established from their non-keplerian rotation curves ([43, 44], 
an innovative review on this issue can be found in [45]). A massive dark component, becoming 
progressively more important at increasing radii and decreasing galaxy luminosity is present in any 
galaxy [46]. This evidence relies on the fact that rotation curves (RCs) of spiral galaxies suitably 
measure their total gravitational potential [see e.g. 47, 48]. More precisely, high-quality RCs measure 
the spiral's circular velocity V{r), i.e. the rotational equilibrium velocity implied by their underlying 
mass distribution. In each galaxy, the former is related to the total gravitational potential 4>tot, which 
is the sum of four distinct mass components: a spherical stellar bulge {B), a Dark Matter halo {DM), 
a stellar disk (*) and a gaseous disc {HI). Such relation follows the equation: 

V\r) = V| + Vl^M + + = rj-^^tot • (3.1) 

where (/)tot = 4^B + 4'dm + 4>* + (j^Hi- It is therefore possible to equal the observed circular velocity, 
i.e the LHS of Eq. 3.1, with the velocity model given by the RHS of Eq. 3.1, thereby inferring the 
radial profile of the total gravitational potential. More specifically, this is obtained by means of the 
Poisson equation from the surface/spatial densities of the various mass components. 

3.1 The Baryonic contribution 

The surface density of the stellar disk S^,(r) is obtained from the observed surface brightness, once 
we assume a radially constant stellar mass-to-light ratio [see 49]. It is generally found [see e.g. 50] 
that this quantity follows the relation: 



where Mr, is the disk mass and i^/p is the disk scale length; therefore, for the RC we get a contribu- 
tion: 

« - (f ) 
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from the stellar disk, where x = r jRo, G is the gravitational constant, and B = IqKq — IiKi is a 
combination of Bessel functions. 

The contribution of the gaseous disk to the RC, V^j, can be directly derived from the HI surface 
brightness; however, for the high-luminosity spirals we consider here, this contribution is always neg- 
ligible [see e.g. 44]. Moreover, although spirals can have a non-negligible bulge velocity component, 
this results virtually indistinguishable from that of the exponential disk [see 51] for the determination 
of the DM structural mass parameters via RC fitting. 

3.2 The phenomenology of RCs of Spirals 

The systematical study of the spiral RCs has highlighted their universal character, with the galaxy 
magnitude M/ playing the role of an additional characteristic parameter identifying specific galaxies 
along a universal behavior of the RC [52, 53]. This has led to the concept of "Universal Rotation 
Curve" (URC), i.e. a function Vjjb.c{t / Rd, Mj) of radius (in units of the disk length scale, r/Ro) 
and of galaxy magnitude (Mj) that well reproduces the RC of any object of known Mj and Rd [see 
52]. In detail, [52] assigned each of their 900 individual RCs of spirals extended out to Ropt — SRd 
to the coiTcsponding magnitude bin among the 11 bins in which the whole spiral /-band magni- 
tude range Mj G [—24.3, —16.3] was divided. The RCs of the objects in each luminosity bin were 
finally averaged in radial bins of width O.lRopt- This led to a family of synthetic coadded RCs 
ycoadd{f / Rd, Ml) [see e.g. Fig. 1 of 52]. A similar result was obtained by [54]. These coadded 
RCs result regular, smooth and with a very small intrinsic variance. Additional data, including very 
extended RCs and the virial velocities Vmr = (GMmr/Rmr)^^'^ [see 55], allowed to extend the coad- 
ded RC out to the galaxies virial radii [53, 55]. With such a general prescription, all kinematical data 
were fit by the same analytical function Vjjrc{t / Rd, Mi) [see again 53]. Further evidence that the 
URC does represent the general RC of any spiral of any luminosity was provided by [54]. 

In this paper we are concerned with RCs of high-luminosity spirals, as their DM halos cor- 
respond to the smallest objects for which the simulations by [34] provide a reliable estimation of 
the halo structural parameters. As a matter of fact, we aim to compare the URC relative to the top- 
luminosity spiral objects - i.e. to objects with M/ < — 23 [see the two last panels in Fig. 1 of 52] 
- with the prediction of some specific mass models. The former is plotted in Fig. 2 as black filled 
circles, with error bars that include also the variance detected object by object in the 56 RCs of [52] 
and the 30 of [56], i.e. of the available samples of high luminosity spiral RCs. More specifically, the 
data points in Fig. 2 show a typical circular velocity profile of a "very big" spiral, M/ = —23.0 and 
consequently with R^ = 7.3 kpc, V{Ropt) = 235km/s [see 53, although no result in the present 
paper depends on these numbers and exactly the same result would arise if we had used anyone of the 
available individual RC of high luminosity spira as all of them are, for the present aim, indistinguish- 
able from the adopted one]. Obviously, any successful mass model must be capable to reproduce 
these data. 

3.3 The velocity profile in ACDM halos 

It is well known that numerical simulations performed in the ACDM scenario lead to virialized halos 
of Dark Matter with a very specific density and velocity profile, called the NFW profile [9] : 

PNFw{l^) = — 2 ■ ^^-^^ 

(r/rs) (1 + r/rs) 

In Eq. (3.4) is a characteristic scale radius and ps is the corresponding density. We remind that 
for the virial radius R^ir and halo mass Myir and the mean universal density pu we have: Myir ~ 
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Figure 2. Circular rotation velocity for the case of a standard ACDM concentration v^. Mass relation (black 
solid line) and for the two cDE models under investigation, EXP010a2 (red solid line) and EXP015a3 (green 
solid line), compared to the URC at Mj = —23 (black points with error bars). 

WOpuR^ir- Furthermore, and Rmr are found to be related within a reasonable scatter according to 
the equation [57] : 

c = R..r/rs ^ 9.7 (^^) . (3.5) 

Although other mildly different density profiles and/or c vs. M^r relations have been proposed, 
adopting Eq. (3.5) as the standard ACDM concentration does not affect the outcomes of the present 
work. 

Since they emerged in N-body simulations, the cuspy NFW density profiles have been in ten- 
sion with those detected around dwarf spirals [16], giving rise to the so called "cusp-core" problem. 
Later, e.g. in [58] and [59], a suitable investigation of a number of proper test-cases was performed 
by means of a careful analysis of 2D, high-quality, extended, regular, and symmetric RCs for spiral 
galaxies. It is well known that in these objects the NFW halo predictions and kinematical data ai^e 
found in strong disagreement in several aspects. In fact, the stellar disk -i- NFW halo mass model 
(NFWD) i) fits very poorly a significant number of RCs; ii) often requires an implausibly low stellar 
mass-to-light ratio and/or ///) an unphysical high halo mass and/or iv) an implausible value for the 
concentration parameter c [see e.g. 59-65]. Up to now, the failure of the NFW velocity profiles is 
very evident in low luminosity galaxies, as their steep RCs cannot be accounted for by a NFWD ve- 
locity profile. Noticeably, instead, for the most luminous objects as large spiral galaxies, their flatfish 
RCs can be still (marginally) reproduced also by a NFWD mass model. As an example, the RCs of 
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Figure 3. Values of x^od^^ best-fit mass models for different assumptions on the concentration parameter 
c/cA. Standard NFW models correspond to c/ca = 1. The Xrcd values for the two cDE models discussed in 
the present work are shown by a red (EXP010a2) and a green (EXP015a3) dot. The latter model provides a 
very poor fit to the data, as shown by the large x^ed value. 

the luminous spirals MW and M31 can be well fitted also by a NFWD mass model [66]. 

The ACDM scenario might generically account for the "cusp-core" tension described above if 
some physical process, related to the galaxy formation process itself, erased the previously formed 
cosmological cusp. A rather common approach is therefore to assume a ACDM scenario and to 
explain the observed RCs postulating that some process has transformed the originally sharp halo 
density profiles into the observed shallow ones. Consequently, any process or scenario in which 
ACDM halos get a more pronounced cusp with respect to the standard case should be disregai^ded. In 
fact, given the present possible tension between ACDM halos predictions and actual observations, it 
is possible to rule out or tightly constrain any scenario that further worsens the discrepancy. In other 
words, any reasonable process beyond the adiabatic formation of CDM halos in ACDM we want to 
invoke must erase the central cusp and not make it more pronounced. 

4 Constraining cDE models with RCs of High-Luminosity Spirals 

In Fig. 2 we plot (as black points) the synthetic RC of spiral galaxies of Mj ~ — 23 from 1 kpc, 
their typical innermost radius with reliable individual measurements, to 330 kpc (their typical virial 
radius). It is worth to recall here that the cored halo -i- Freeman Disk model (not discussed here) 
perfectly fits these data by passing all the data points within their error bars [53]. 

It is also not surprising that the NFWD model with the standard concentration relation (3.5) fits 
the same data in a quite satisfactory way (black line in Fig. 2 and black filled circle in Fig. 4 where 
the resulting reduced chi-squared Xred shown). In fact, so far in the literature there are no cases in 
which the flattish individual RCs of high-luminosity spiral galaxies can be unambiguously modeled 
exclusively by a cored DM distribution as it occurs for the steeper RCs of lower luminosity objects. 

Let us now investigate the case - as for the cDE scenarios with a steep growth in time of the 
coupling function Pc [34] - in which the emerging dark matter halos have a NFW density profile 
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with a value for the concentration parameter significantly higher with respect to what prescribed by 
Eq. (3.5) for halos of this mass For instance, we can consider concentration values as found in N-body 
simulations for the EXP010a2 and EXP015a3 models investigated in the present work (see Fig. ??), 
which at masses of the order of a few x 10^^ Af© are given by: 

EXP010a2 CW1.3-CA, (4.1) 
EXP015a3 cRil.85-CA, (4.2) 

where c\ is the standard value of the concentration predicted by ACDM simulations and consistent 
with Eq. 3.5. 

With these values of the concentration parameter, we can now attempt fitting the reference RCs 
of high-luminosity spirals shown in Fig. 2 with a NFWD mass model for NEW DM halos of different 
concentration, and test the quality of the fit with the reduced chi-squared Xrcd defined as: 

2 = 1 V^URC - K.od 4 3 

Xrcd - ^ _ 2 0.0175 • V^j^c ' 

where N=19, the quantity in the denominator is the observational enw of VuRC and Fmod is the 
velocity model under investigation. For a standard ACDM model, coiTcsponding to c/ca = 1, the 
best fit parameters are Mmj. = 1.1 x 10^^ Mq, Md = 5.1 x IO^^Mq with a reduced chi-squared of 

The fit appears instead clearly unsuccessful (see the green line in Fig. 2 and the green point 
in Fig. 4) for the most extreme cDE case EXP015a3, corresponding to the parameters (/3o,/3i) = 
(0.75,3). The best fit values of the parameters for this model are Mi,ir = 6 x 1O^^M0, Md = 
2.6 X 1O^^M0, but the fit results extremely poor with a reduced chi-squared of Xrcd ~ 
the other hand, the fit is not significantly worse than in the standard ACDM case for the other, less 
extreme cDE scenario discussed in this work, EXP010a2, with (/3o, /3i) = (0.5, 2), for which the best 
fit parameters are M^i,. = 8.7 x IO^^Mq, Md = 4.1 x lO^U/© with a x^cd ~ l-'^- Let us stress 
here that the fitting uncertainty in the the two structural parameters is large but not relevant here. In 
the former case, no model can fit the data in a satisfactory way. 

More in general, let us consider DM halos with a NEW profile but with a c-M relation different 
from the standard ACDM case, i.e. c/ca ranging from c/ca = 0.8, leading to a mai^ginally shal- 
lower density profile, to c/c\ = 2, i.e. a significantly more concentrated DM halo. We attempt in 
these cases a NFWD modeling of the reference RC. The result is clear (see Fig. 4): at a value of 
c/cA ~ 1.8 the fit becomes very problematic and even flat RCs rule out a NEW halo with such a high 
concentration. 



5 Discussion and Conclusions 



In the ACDM framework, DM virialized halos are characterized by a cuspy density distribution that 
is certainly not seen in the halos around dwarf spiral galaxies and that might be absent also in high- 
luminosity spiral galaxies. Some alternative cosmological models, motivated by the need to address 
the fundamental problems of the cosmological constant, as e.g. some specific models of interaction 
between Dark Energy and Dark Matter can modify the typical concentration of DM halos as a function 
f halo mass, both reducing it or increasing it depending on the specific form of the DE self-interaction 
potential and of the coupling function. However, while the reduction of halo concentration, when it 
occurs, is in general quite modest for presently viable cDE models, some specific scenarios can 
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lead to DM halos that are significantly more concentrated than the standard ACDM ones. As a 
consequence of this, also the rotation curves of high-luminosity spiral galaxies as obtained from the 
Universal Rotation Curve or from the kinematics of numerous suitable objects, cannot be reproduced 
in a satisfactory way by a mass model that accounts for the effect of enhanced concentration in 
interacting DE cosmologies. 

In the present work, we have investigated the impact that these specific realizations of the in- 
teracting DE scenario have on the rotation curves of luminous spiral galaxies. More specifically, we 
have restricted our analysis to a particular- class of coupled DE scenarios characterized by a steep 
growth in time of the interaction strength, parametrized as a positive power of the cosmic scale fac- 
tor. Such models have been shown to be in full agreement with present constraints on the background 
evolution of the Universe, as they do not feature the typical early DE scaling of other coupled DE sce- 
narios like e.g. models with a constant or a more slowly evolving coupling function. The background 
evolution of the models investigated in this work is in fact practically indistinguishable from the stan- 
dard ACDM cosmology. Furthermore, the impact of these alternative cosmologies on the statistical 
properties of large scale structures is confined at very low redshifts, which makes it problematic to 
test the models via e.g. weak lensing measurements. 

However, as mentioned above, a series of high-resolution hydrodynamical N-body simulations 
for this specific class of interacting DE cosmologies have recently shown that the fast growth of 
the coupling strength determines a significant increase of the normalization of the concentration v^^. 
Mass relation for CDM halos at low redshifts, thereby allowing for a direct test of the models through 
dynamical probes in collapsed structures. In the present study we have therefore considered the 
effect that an average higher halo concentration at halo masses in the range of luminous spirals - 
consistent with the results of N-body simulations for two specific realizations of coupled DE models 
- determines on the rotation curve of the galaxy, and compared the inferred circular velocity with 
real data, namely with the Universal Rotation Curve that for our purpose well represents the typical 
circular velocity of spirals. 

While the standard ACDM value of the halo concentration provides a marginally acceptable 
fit to the data, the increase of the concentration due to the interaction between DE and DM for the 
specific models investigated in this work determines a significantly worse fit and exacerbates the 
problem of the central Dark Matter cusps in galactic halos. In particular, our results show that for 
the most extreme of our coupled DE scenarios the best fit mass modeling has a reduced chi-squared 
of Xred ~ thereby allowing to rule out the model. The other model under investigation, which 
features a smaller present value and a shallower growth in time of the coupling function, shows 
instead a still marginally acceptable fit to the data. 

We can therefore conclude that the nonlinear- effects of an interaction between DE and CDM 
provide a direct way to constrain the evolution of the coupling function, and that the parameter space 
of coupled DE models can be constrained by means of direct observations of the dynamical properties 
of galaxies. Furthermore, since the problem of the central cuspiness of DM halos is mostly severe 
for faint and low-mass objects, we argue that a detailed investigation of the effects of coupled DE 
models on the halo concentration at a lower mass range than allowed by the present resolution of our 
simulations might result in tighter constraints on the coupled DE parameter space. It is also important 
to stress here that a large fraction of the parameter space of coupled DE cosmologies does not result 
at all in an enhancement of halo concentrations, with a large number of models even showing a slight 
decrease of the concentrations over a wide range of masses. These realizations of the coupled DE 
scenario, not discussed in the present work, are clearly not affected by our constraints and should 
still be considered as viable alternatives to the ACDM paradigm. Nevertheless, our results show for 
the first time that it is possible to constrain the properties of Dark Energy also through the dynamical 
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features of luminous galaxies. 
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